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Abstract: 

'--^The  Influence  of  strong  multiple  interplanar  forces  is  shown  to  give 
rise  to  a  new  class  of  surface  vibrational  resonances  on  metal  surfaces  as 
a  result  of  higher  harmonic  terms  In  the  bulk  phonon  dispersion  relations. 
This  effect  is  Illustrated  with  the  bcc  Fe(lll)  surface  where  a  surface 
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With  recent  advances  In  the  experimental  techniques  for  the  study  of 
surface  vibrations,1  considerable  Interest  exists  In  determining  the 
relation  between  the  geometric  structure  of  the  crystal  surface  and  the 
accompanying  surface  vibrational  modes.2  Recent  electron  energy  loss 
measurements  have  shown  the  existence  of  a  new  type  of  surface  vibrational 
resonance  on  the  fcc(llO)  surfaces,3  which  resulted  from  the  unique 
coordination  of  the  atoms  along  the  fcc(llO)  direction.  This  particular 
coordination  along  the  fcc(llO)  direction  results  In  a  nonmonotonic 
longitudinal  bulk  phonon  dispersion  and  a  pseudo  band  gap  which  give  rise 
to  an  eigenmode  localized  at  the  surface  forming  a  vibrational  resonance. 

In  this  paper,  we  show  that  this  new  class4  of  surface  vibrational 
resonances  Is  a  general  effect  which  results  from  higher  harmonic  terms  In 
the  bulk  phonon  dispersion  relation  for  the  longitudinal  phonons 
propagating  normal  to  the  surface.  These  higher  harmonic  terras  allow  the 
existence  of  a  localized  solution  at  the  surface- vacuum  Interface  which 
gives  rise  to  a  surface  vibrational  resonance.  We  Illustrate  the 
generality  of  this  new  class  of  surface  resonances  with  a  lattice  dynamical 
analysis  of  the  Iron  bcc ( 111 )  surface  where  this  type  of  surface  resonance 
Is  observed  by  electron  energy  loss  spectrosocpy  (EELS).  The  case 
considered  here  corresponds  to  the  r  point  of  the  surface  Brlllouln  zone  In 
order  to  compare  to  the  Inelastic  dipole  scattering  observed  In  the  EELS 
measurements^ 6  Our  calculations  show  that  these  resonances  disperse  across 
the  surface  Brillouin  zone  and  should  also  be  observed  at  higher  parallel 
wavevector  transfers.7 

The  EELS  measurements  were  performed  In  a  multi  technique  ultrahlgh 
vacuum  system  with  a  base  pressure  of  4* 10" 11  torr,  which  Is  described  In 
detail  elsewhere.8  The  EELS  spectrometer  consists  of  a  double  pass  127* 


3 


cylindrical  deflector  for  both  the  monochromator  and  analyzer.  The 
electron  Impact  energy  was  2.9  eV.  The  scattering  plane  containing  the 
Incident  and  scattered  electrons  was  defined  by  the  crystal  normal  and  the 
crystal  [Oil]  direction.  The  crystal  was  cleaned  extensively  In  previous 
studies.9  The  surface  cleaning  procedure  for  the  present  study  Included 
Neon  ion  sputtering  at  1  keV  and  annealing  to  850  K.  Surface  cleanliness 
was  monitored  with  Auger  electron  spectroscopy  and  EELS. 

The  longitudinal  phonon  (dynamics  along  high  symmetry  directions  In 
cubic  crystals  reduces  to  a  well-known  one-dimensional  problem  as  a  result 
of  the  factorization  of  the  dynamical  matrix.  The  one-dimensional  dynamics 


Involves  the  rigid  displacement  of  lattice  planes,  where  the  equations  of 
motion  are  given  by,10*11 


0)2U  =  -l  —  (u  -  u  ) 


where  u*  Is  the  displacement  of  an  atom  In  layer  t  normal  to  that  layer, 

M  Is  the  mass  of  the  atom,  <■>  Is  the  phonon  frequency,  and  +n  Is  the  force 

on  a  given  atom  from  all  the  atoms  In  the  neighboring  plane  n.  For  the 

Infinite  lattice  the  equations  of  motion  admit  simple  plane  wave  solutions 

where  the  resulting  phonon  dispersion  is  given  In  terms  of  a  Fourier  series 

In  the  Interplanar  force  constants  $n  as, 

2$ 

0>2U)  »  l  (l  -  cosnwc),  (2) 

n  M 


where  c  Is  the  reduced  wavevector,  c  *  Qzd/w,  and  d  is  the  Interplanar 
distance.  This  relation  can  then  be  used  to  determine  the  Interplanar 
force  constants  from  the  measured  bulk  phonon  dispersions  by  Inelastic 
neutron  scattering  as  has  been  done  extensively  In  the  past.12 


Surface  modes  occur  In  the  semi- Infinite  crystal  as  a  natural 
consequence  of  the  broken  translational  symmetry  normal  to  the  surface. 

For  the  semi-infinite  crystal  the  solutions  to  the  equations  of  motion  In 
Eq.  (1)  are  no  longer  simple  plane  waves  but  are  now  modified  by  the 
reflection  of  the  phonons  at  the  surface.3  The  surface  can  be  viewed  as  a 
perturbation  which,  given  an  Incident  wave  on  the  surface,  will  mix 
together  reflected  waves  with  the  same  energy  and  parallel  wavevector.  A 
special  situation  occurs  when  the  series  In  Eq.  (2)  contains  terms  with 
n  ^  2.  For  n  >.  2  the  wavevector  corresponding  to  the  reflected  phonons 
becomes  complex  In  certain  frequency  ranges  as  determined  from  the 
dispersion  relations  In  Eq.  (2).  These  complex  wavevectors  then  gives  rise 
to  localized  evanescent  phonons  at  the  surface.  The  situation  where  n  *  2 
was  observed  on  the  fcc(llO)  surface.3  However,  on  the  bcc(lll)  surface 
the  interplanar  interactions  are  more  extensive  and  terms  at  least  up  to  n 
*  3  need  to  be  included  In  order  to  describe  the  lattice  dynamics. 

The  bulk  longitudinal  phonon  dispersion  along  the  Fe(lll)  direction  Is 
shown  In  Fig.  1,  as  measured  by  inelastic  neutron  scattering.13  A 
characteristic  feature  of  the  large  higher  harmonic  terms  In  the  phonon 
dispersion  Is  the  appearance  of  additional  maxima  and  minima  within  the 
Brillouin  zone  as  seen  In  Fig.  1.  The  occurrence  of  strong  multiple 
interlayer  forces  for  the  bcc(lll)  direction  may  be  seen  from  examining  the 
atomic  structure  along  this  direction.  As  Fig.  2  shows,  an  atom  In  a  layer 
parallel  to  the  bcc(lll)  surface  has  nearest  neighboring  atoms  In  the  first 
and  third  neighboring  layers.  Next  nearest  neighboring  atoms  are  located 
In  the  second  neighboring  layer.  Thus  Including  only  nearest  and  next 
nearest  neighboring  forces,  the  dispersion  along  this  direction  will  have 
up  to  three  harmonic  terms  in  the  series  expansion  In  Eq.  2  given 


explicitly  as, 

w2(?)  a  1  {4»i(l  -  COSIT?)  +  $2^1  "  COS2ir?)  +  ^(l  -  C0S3ir?)}.  (3) 

M 

The  interlayer  force  constants  may  be  obtained  from  a  previous  general 
Born-von  Karman  analysis  of  the  neutron  data,13  and  including  only  the 
nearest  and  next  nearest  interactions  they  are  given  by:  ♦1/M  *  96.38, 

♦2/M  =  68.35,  and  ^/M  *  219.14  meV2.  The  dispersion  obtained  with  these 
force  constants  is  in  good  agreement  with  the  neutron  data  as  shown  by  the 
solid  line  in  Fig.  1.  From  the  bulk  dispersion  we  see  that  there  are  three 
propagating  solutions  for  bulk  phonons  having  energies  in  the  range  23  <  hu 
<  35  meV.  For  energies  <  23  meV  two  of  the  solutions  will  exhTBTt 'complex 
wavevectors.  For  the  semi -infinite  crystal  these  complex  solutions 
correspond  to  evanescent  phonons  localized  at  the  surface.  Since  we  still 
have  a  single  propagating  solution  for  this  range  of  frequencies,  these 
evanescent  phonons  can  couple  to  the  overlapping  bulk  phonons  and  form  a 
surface  resonance.  This  effect  may  be  seen  explicitly  by  examining  the 
phonon  density  of  states  for  the  semi-infinite  crystal  projected  on  the 
motion  of  the  outer  surface  layers. 

The  phonon  density  of  states  for  the  semi-infinite  crystal  may  be 
evaluated  by  using  an  efficient  transfer  matrix  technique.14  Figure  3 
shows  the  phonon  density  of  states  for  the  longitudinal  phonons  in  Fig.  1 
projected  on  the  motion  of  a  bulk  layer  (upper  panel)  and  the  outer  surface 
layer  (lower  panel).  The  density  of  states  for  the  bulk  layer  shows  a 
region  of  high  density  of  states  for  energies  >23  meV,  with 
multiple  van  Hove  singularities  due  to  the  three  stationary  points  in  the 
phonon  dispersion  in  Fig.  1.  For  energies  <  23  meV  the  bulk  states  are 
severely  depleted  and  hence  this  region  is  called  a  pseudo  band  gap.3*7 
In  this  projection  the  bulk  density  of  states  tends  to  a  constant  value  at 


low  energies,  which  Is  simply  proportional  to  the  sound  velocity  In  the 
dispersion  shown  In  Fig.  1. 


The  lower  panel  In  Fig.  3  shows  the  density  of  states  projected  on  the 
motion  of  the  outer  surface  layer.  The  surface  force  constants  were 
obtained  from  considering  the  loss  In  coordination  obtained  from  a 
truncation  of  the  bulk  crystal.  The  interlayer  force  constants  were  then 
simply  taken  from  the  bulk  values  determined  above.  As  seen  in  this 
surface  projection  a  very  sharp  peak  Is  observed  within  the  bulk  pseudo 
band  gap  at  19.7  meV.  Conceptually,  the  formation  of  this  resonance  Is 
similar  to  the  case  of  true  surface  phonons  which  are  observed  to  be  split 
off  from  the  bottom  of  bulk  manifolds  Into  an  absolute  bulk  band  gap.5  In 
the  present  case  we  have  a  pseudo  band  gap  with  a  depleted  but  finite 
density  of  bulk  states.  The  surface  then  splits  off  a  mode  from  the  region 
of  high  density  of  bulk  states  and  becomes  resonant  with  the  bulk  states  In 
the  pseudo  band  gap.  In  comparison  to  the  fcc(llO)  surfaces,  the  resonance 
for  the  bcc(lll)  surface  Is  significantly  stronger  and  sharper.  The 
Intrinsic  width  for  the  Fe(lll)  surface  Is  0.9  meV  compared  to  3.5  meV 
calculated  for  the  NI(llO)  surface.3  This  results  from  the  enhanced 
depletion  of  bulk  states  In  the  pseudo  band  for  the  bcc(lll)  structure  In 
comparison  to  the  fee (110)  structure.  The  additional  small  structure  at 
21.3  meV  In  the  surface  density  of  states  is  due  to  a  second  resonance 
which  is  mainly  localized  on  the  second  and  third  layers,  while  the  former 
at  19.7  meV  Is  mainly  localized  on  the  first  and  second  layers. 

A  comparison  with  the  predictions  of  the  surface  lattice  dynamical 
calculations  can  be  made  with  the  measured  EELS  loss  spectra  shown  In  Fig. 
4.  A  loss  feature  attributed  to  the  structural  surface  resonance  is 
measured  at  21  meV.  The  small  structure  at  -  60  meV  Is  due  to  a  small 
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amount,  <  1%,  of  oxygen  contamination  which  does  not  affect  the  observation 
of  the  surface  resonance  for  very  low  oxygen  coverages.  The  measured  EELS 
Intensity  for  the  dipole  active  losses  Is  directly  proportional  to  the 
surface  density  of  states,  g(u),  projected  onto  a  normalized  dipole  active 
charge  field  v*,14 

g{w,{v*})  »  l  |£  v*5^|2  <5(u  -  o>x),  (4) 

where  Is  the  normalized  displacement  of  the  substrate  layer  t  for  an 

«• 

eigenmode  x  with  energy  Ifox  of  the  dynamical  matrix  D.  For  the  case  of  a 
clean  metal  surface  considered  here,  the  dipole  projection  Is  given  by  the 
projection  on  the  relative  motion  of  the  outer  two  surface  layers.3*15 
The  dipole  projected  density  of  states  Is  shown  in  the  Inset  In  Fig.  4  for 
a  comparison  with  the  measured  loss  spectra.  The  density  of  states  was 
convoluted  with  a  4.5  meV  Gaussian  distribution  to  simulate  the 
instrumental  linewidth.  The  dipole  projection  gives  a  spectral  loss  peak 
at  20.6  meV  which  is  in  excellent  agreement  with  the  measured  loss 
spectra.  The  slight  shift  In  the  dipole  projection  from  the  19.7  meV 
calculated  for  the  outer  surface  layer  projection  Is  due  to  the  suppression 
of  the  asymmetry  in  the  Fano-like  shape  of  the  resonance  (see  Fig.  3)  In 
the  dipole  projection.  In  the  dipole  projection  shown  in  Fig.  4,  the 
density  of  states  goes  to  zero  at  low  energies,  since  the  two  outer  surface 
layers  are  moving  in  phase  In  the  long  wavelength  limit. 

In  conclusion,  we  have  shown  that  strong  multiple  Interlayer  forces 
give  rise  to  higher  harmonic  terms  in  the  bulk  phonon  dispersion  relations 
which  lead  to  strongly  depleted  pseudo  band  gaps  and  evanescent  phonons 
localized  at  the  surface  producing  vibrational  resonances.  This  type  of 
resonance  is  most  noticeable  In  crystal  structures  which  have  nearest 


neighboring  forces  In  multiple  neighboring  planes,  as  Is  Illustrated  by  the 
bcc(lll)  structure. 
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FIGURE  CAPTIOUS 

Fig.  1.  The  iron  [111]  bulk  longitudinal  phonon  dispersion.  (  ■  ): 

Measurements  by  inelastic  neutron  scattering.13  ( - ):  The 

calculated  dispersion  using  the  first  three  interlayer  force 
constants,  including  only  nearest  and  next  nearest  neighboring 
forces,  obtained  from  a  previous  general  Born- von  Karman  analysis 
of  the  neutron  data.13  The  force  constants  are:  4^/M  =  96.38, 
♦2/M  *  68.35,  and  *3/M  =  219.14  meV2. 

Fig.  2.  Schematic  diagram  of  the  bcc(lll)  surface  structure,  (a)  Top 

view.  The  numbers  denote  sequential  layers  of  the  substrate  with 
number  one  corresponding  to  the  surface  layer,  (b)  Side  view. 
Notice  an  atom  has  nearest  neighbors  in  the  first  and  third 
neighboring  layers.  Next  nearest  neighboring  atoms  are  in  the 
second  neighboring  layer. 

Fig.  3.  Vibrational  density  of  states  for  the  longitudinal  phonons 

propagating  normal  to  the  surface.  The  upper  panel  shows  the 
projection  on  the  motion  of  a  bulk  layer.  The  lower  panel  shows 
the  projection  on  the  motion  of  the  outer  surface  layer.  Two 
resonances  are  observed  in  the  surface  density  of  states  at  19.7 
and  21.3  meV.  The  resonance  at  19.7  meV  is  mainly  localized  on 
the  first  and  second  layers,  while  the  latter  at  21.3  meV  is 
mainly  localized  on  the  second  and  third  layers. 

Fig.  4.  An  electron  energy  loss  spectra  of  the  Fe(lll)  surface.  The 
spectrum  was  recorded  at  300  K  in  the  specular  direction.  The 
surface  resonance  is  observed  In  the  loss  spectrum  at  21  meV.  The 
small  structure  at  "  60  meV  Is  due  to  a  small  amount,  <  1%,  of 
oxygen  contamination.  The  inset  shows  the  dipole  projected 
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density  of  states  for  a  comparison  with  the  measured  loss 
spectra.  The  projection  Is  taken  on  the  relative  motion  of  the 
outer  two  surface  layers  and  Is  convoluted  with  a  4.5  meV  Gaussian 
distribution  to  account  for  the  Instrumental  resolution. 
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